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Abstract 
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and to excitation energies of 2^ states in '^^Ca and ^^S mirror nuclei are studied using the real- 
energy continuum shell model. 
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I. INTRODUCTION 



Atomic nuclei form a network of coupled systems communicating with each other through 



decays and captures [l|. If continuum states are neglected, this communication is broken and 
each system becomes an isolated closed quantum system (CQS). It is obvious, that the CQS 
description of atomic nuclei (e.g. the nuclear shell model (SM)) becomes self-contradictory 
for weakly-bound or unbound states. 



A classic example of a continuum coupling is the Thomas- Ehrmann shift {2] which man- 
ifests itself in the asymmetry in the energy spectra between mirror nuclei having different 
particle emission thresholds. A consistent description of the interplay between scattering 
states, resonances, and bound state requires an open quantum system (OQS) formulation. 
Comprehensive many-body theory of weakly bound/unbound states has been advanced re- 
cently in the time-asymmetric quantum mechanics using the complete ensemble of single- 
particle states consisting of resonant (Gamow) states and the complex-energy, non-resonant 
continuum of scattering states from which the complete many-body basis of OQSs can be 
obtained Another formulation of the continuum shell model is obtained by embed- 
ding standard SM in the continuum of decay channels. This approach provides a unified 



description of nuclear structure and nuclear reaction aspects 

In this paper, we study the effect of the continuum coupling on spectra of ^^Ca and 
^^S mirror nuclei. We show that the continuum coupling explains naturally not only the 
appearance of the asymmetry in spectra but also provides a major part of the Thomas- 
Ehrmann shift. We shall also discuss salient effects of a continuum coupling in the binding 
energy systematics, in particular the anti-odd-even staggering (anti-OES) and the effective 
range of energies in which various decay channels are correlated with each other with the 
discrete many-body states changing significantly their energy and wave function. 

II. SHELL MODEL EMBEDDED IN THE CONTINUUM 

^ In the shell model embedded in the continuum (SMEC), nucleus is described as an OQS 
. The total function space consists of the set of £^-functions, as in the standard SM, 
and the set of scattering states (decay channels). These two sets are obtained by solving 
the Schrodinger equation for discrete states of the closed subsystem (closed quantum sys- 
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tern (CQS)): HsM^i = E\^^^^i , and for scattering states of the external environment: 
— Hcc')Ce^^ — ; where Hsm is the SM Hamiltonian, and H^^ = Hq + V^^ is the 
coupled-channel (CC) Hamiltonian. Channels: c = [J/^~^; (Z, j)]"'**, are determined by the 
motion of an unbound particle with orbital angular momentum I and total angular momen- 
tum j relative to the residual nucleus with A — \ bound particles in a SM state <l>/-\ 
are channel projected scattering states with outgoing asymptotics. States of the daughter 
nucleus are assumed to be stable with respect to the particle emission. By means of two 
functions sets: Q = }, V = {C^;^"*}) can define the corresponding projection op- 
erators and the projected Hamiltonians: QHQ = Hqq and PHP = Hpp jj, [sj. Hsm is 
identified with the CQS Hamiltonian Hqq and H^c with Hpp. The coupling term HpQ is 
given by the two-body residual interaction |8[. 

Schrodinger equation in the total function space splits into the two equations for projected 
operators Hqq, Hpp. Assuming Q + V = Id, one can determine a third wave function uj^^^^ 
which is a continuation of SM wave function $^ in the scattering continuum, uj^^^ is obtained 
by solving the CC equations with the source term j^. Using the three function sets: }, 
{Ce^^} and {co'l^'*}, one obtains a solution in the total function space: 

n = CE + Y.^^t + ujt\Emt\{E - nQQ{E))-'\^t){^k\HQP\CE) (i) 

where E is the total energy and 1-Lqq{E) is the energy- dependent effective Hamiltonian in 
Q subspace: 

nQQ{E) = HQQ + HQpG'^;\E)HpQ , (2) 

where G^j^\e) in ([2]) is the Green function in V. Ti-QQ takes into account a modification of 
the CQS Hamiltonian {Hqq) by couplings to the environment of decay channels. Hqq is a 
complex-symmetric matrix above the particle-emission threshold E^'^^'^^ and Hermitian below 
it. Diagonalization of ?^qq by an orthogonal and, in general, non-unitary transformation 
yields complex eigenvalues Ei — ^iTi, which depend on the energy E of the particle in the 
continuum. Energies and widths of the resonance states follow from: E^ = Ei{E = E^), 
Fj = Ti{E = Ei), where Ei{E) and ri(-E) are the eigenvalues of T-Cqq{E). Details of the 
SMEC calculations can be found in Refs. U, 
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III. FEATURES OF THE CONTINUUM-COUPLING ENERGY CORRECTION 
TO EIGENVALUES OF THE CLOSED QUANTUM SYSTEM 

In this chapter, we shall discuss salient features of the continuum-coupling energy cor- 
rection: e\^°"\e) = {^i\HQpG^p\E)HpQ\^i) , for the ground state (g.s.) of neutron-rich 
oxygen and fluorine isotopes. Details of the SMEC calculations for continuum-coupling 
energy correction to binding energies can be found in Ref. In these calculations, all 
possible couplings of the g.s. of A nucleus to the states of A — 1 nucleus are taken into 
account incoherently. 



A. Binding systematics in neutron-rich nuclei 

For T = 1 couplings, the average behavior of EgT^^^ is determined by the strong de- 
pendence on En^"^^ (cf Fig. la of Ref. [6] for oxygen isotopes)^. Close to the neutron 
drip line, this leads to an effective enhancement of nn continuum-coupling strength which 
cannot be compensated by the -En '^'^■'-independent correction of monopole terms. On the 
top of this behavior, one can see the odd-even staggering (OES) of E^°''\n). Blocking 
of the virtual scattering to continuum states by an unpaired nucleon diminishes the nn 
continuum-coupling energy correction in odd-A^ nuclei. For a flxed En^'^\ the coupling of 
[A^ = 2k, Z] system to [N = 2k — 1, Z] ^ n decay channels is enhanced and the coupling 
of [N = 2k + 1, Z] system to [A^ = 2k, Z]®n decay channels becomes weaker with respect 
to an averaged behavior. This drip-line effect is seen only in a narrow range of excitation 
energies around E^n'^'\N) ~ and vanishes for E^n''\N) > 4 MeV. For more realistic values 
of separation energies, as given by nuclear SM [6], En^"^^ exhibits the pairing induced OES 
which becomes partially compensated close to a drip line by the anti-OES effect induced by 
couplings to decay channels, both opened and closed. 

In fluorine isotopes (cf Fig. 2a of Ref. Ql), E^°^^^ is dominated by up continuum couplings 
(T = 0, 1). The strength of up continuum coupling can be deduced by comparing SMEC 
results for binding energies with experimental data [6]. An optimal value of V^'^'^^ /V^^"'^ 
close to the neutron drip line is Vq'^^^ ~ {1 /2)Vq^"\ whereas for nuclei close to the valley of 
beta-stability a standard choice is: Vq'^^^ ~ 2Vq"^\ Since the np couplings provide a major 

^ In the chosen sdfp model space, np continuum couphngs are absent in oxygen isotopes. 
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part of this correction, therefore the gradual reduction of Vq"'^^ /Vq'^"'^ toward the neutron 
drip hne leads to a non-linear dependence of the E^°^^^ and, hence, to the dependence of 
two-body monopole terms on the neutron number. A similar dependence is expected if 
effective three-body interactions are included in the two-body framework of the SM 9(]. 

The np continuum coupling in odd-odd {[N = 2k + 1, Z = 2m + 1]) fluorine isotopes is 
increased as compared to the neighboring even-odd {[N = 2k + 2^Z = 2m + 1] and [A^ = 
2k, Z = 2m + 1]) nuclei. Contrary to the continuum-coupling correction for like particles 
{nn or pp couplings), the characteristic anti-OES of E^""^"^^ persists even for M*^'^ > 4 MeV. 
This effect is further enhanced by the OES of one-neutron (In) emission thresholds which 
yields lower En^^^ and, hence, larger continuum-coupling correction for odd-A^ systems. The 
np continuum-coupling energy correction attenuates the OES and can even wash it out close 
to drip lines if the ratio Vq^^^ /Vq"'"'^ would not be strongly reduced from its accepted value 
~ 2 in well-bound nuclei. 

In BCS formalism, the OES is associated with the blocking of a quasi-particle state close 
to the Fermi energy by an unpaired neutron ( resp. proton) what weakens nn (resp. pp) 
pairing correlation in odd-A^ (resp. odd-Z) isotopes. The proximity of continuum states 
in weakly-bound nuclei makes the blocking mechanism less effective, reducing the OES of 
one-nucleon separation energies. This reduction appears even though the strength of nn 
(pp) pairing correlations is essentially unchanged. 



B. Anatomy of the continuum-coupling correction 

A typical behaviour of the total continuum-coupling correction E^°^^^ to the g.s. SM 
energy (the CQS eigenvalue) for different oxygen isotopes is presented in Fig. [T]as a function 
of the neutron energy. This correction is an incoherent sum of contributions from couplings to 
all SM states in the daughter nucleus. E = corresponds to a position of the first In emission 
threshold. Rapid change of Ecorr related to opening of next In emission threshold can be 
seen for ^^O at ~ 3 MeV. The continuum-coupling energy correction rises with number 
of valence neutrons and in general is bigger in even-A^ isotopes. Few notable exceptions 
can be seen right after the closure of the sd shell (cf ^^O in Fig. [1]) and at the begining 
of the sd shell. Behavior of the continuum-coupling correction depends on / of the neutron 
wave involved in a decay channel [Ji'^~^', (^7 j)]'^*^'^' • This correction is largest exactly at 
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FIG. 1: The continuum-coupling energy correction to the SM g.s. energy for different oxygen 
isotopes is plotted as a function of the neutron energy. Contributions from couplings to all available 
states in the daughter nucleus are added incoherently. 

the threshold of a channel [Ji'^^^] i}i3)Y^'''' only for / = neutrons. For higher /-values 
or for protons, the centrifugal barrier and/or the Coulomb barrier shift the maximum of 
energy correction above the threshold. This can be seen for ^^O, where I = 2 neutron wave 
dominates in the g.s. to g.s. coupling. 

Couplings to excited states in A — 1 nucleus are less important in even-iV isotopes. On 
the contrary, in odd-A^ isotopes they dominate. The distribution of contributions to i^gs""^^ 
as a function of the energy of corresponding states in A — 1 nucleus, reflects the nature 
of pairing correlations in odd-A^ and even-A^ isotopes. Strong dissimilarity of g.s. to g.s. 
couplings between odd-A^ and even-A^ oxygen isotopes is attenuated if couplings to decay 
channels involving excited states in A — 1 nucleus are taken into account. The fraction of 
^(corr) (.Qj^jj^g fj-om coupling to the g.s. of daughter nucleus is shown in Figs. [2] and [3] for 
oxygen and fluorine isotopes, respectively. In oxygen isotopes (T = 1 couplings), one can 
see a strong OES following the staggering of neutron pairing correlations in these isotopes. 
In general, strong nn pairing correlations in even-A^ isotopes increase the weight of a g.s. 
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FIG. 2: The fraction of the total continuum-coupling energy correction E^"^^"^ in oxygen isotopes, 
coming from a coupling to the g.s. of a daughter nucleus. SMEC calculations are performed for 
fixed In-threshold energies: En^^^ = 0, 4 MeV. 

to g.s. contribution. On top of this effect, one sees different values of the g.s. fraction at 
the beginning of the shell (filling of 6/5/2 subshell) and at the end of the shell (filling of si/2 
and d3/2 subshells). This picture changes qualitatively in fluorine isotopes (cf Fig. [3]). As 
pointed out in Sect. IIIIl np continuum coupling dominates in this isotopic chain, removing 
most of the OES due to the nn correlations. As in oxygen isotopes, the fraction of E'gs™'^^ 
coming from coupling to the g.s. of A — 1 nucleus is smaller in ^5/2 subshell than in si/2 
and (i3/2 subshells. Maximal value of this fraction is less than ~ 20% in fluorine isotopes, 
whereas for ^^O, ^^O and ^^O the g.s. fraction is > 60%. 



IV. MIRROR-SYMMETRY BREAKING EFFECT OF A CONTINUUM COU- 
PLING: EXAMPLE OF 2+ STATES IN ^^CA AND ^^S 

Symmetry breaking effects in spectra of mirror nuclei are directly or indirectly related to 
the Coulomb interaction. Direct effects (the Thomas- Ehrmann effect) have been extensively 
discussed in the literature. Much less known indirect effects of the Coulomb interaction 
result from different positions of n/p thresholds which modify continuum couplings effects 
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FIG. 3: The same as in Fig. [2] but for fluorine isotopes. 



and, hence, change spectra in mirror systems 8|]. Recently, the energy of 2^ state has been 
measured in a weakly-bound ^^Ca [10]. An excitation energy of this proton-unbound state 
differs by ~ 240 keV from an energy of the well-bound mirror state in ^^S. In this section, 
we will discuss in SMEC the mirror-symmetry breaking continuum-coupling effects for these 
states. Details of SMEC calculations and the choice of an effective interaction are the same 
as described in Sect. Illlland in Ref. {g]. 

SMEC excitation energies of 2^ states in '^^Ca and ^^S are presented in the third column 
of Table [11 The difference of excitation energies in those mirror states is close to the experi- 
mental value (cf the first column of Table [11), but absolute excitation energies are ~ 0.5 MeV 
bigger than experimental ones and ~ 1 MeV bigger than a SM value (cf the second column 
of Table [T]). One should notice that proton continuum couplings for '^^Ca and neutron con- 
tinuum couplings for ^^S result in a small difference of 2f excitation energies in ^^Ca and 
^^S. Almost a whole difference of 2^ excitation energies is due to the neutron-continuum 
couplings for ^^Ca and proton-continuum couplings for ^^S. This is due to the same struc- 
ture of protons in ^^Ca and neutrons in ^^S for both 0^ g.s. and 2^ excited states (the shell 
closure). 

Let us take a closer look at the individual contributions from couplings via continuum to 
different states in A — 1 nuclei. In Fig. [H contributions to the continuum-coupling energy 
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TABLE I: Excitation energy of 2^ states in ^°Ca and S nuclei. Column labelled 'SM' contains 
the SM results. The following three columns present SMEC results with different couplings. The 
third column shows results with both neutron and proton continua included. The next columns 
give the 2]*" energy with only proton or neutron continuum couplings included. All energies are in 
MeV. 



Exp(36Ca) 


SM(36Ca) 


SMEC(36ca) 


36Cap5K + p] 


36Cap5ca + n] 


3.05 ± 0.05 


2.641 


3.403 


3.019 


3.514 


Exp('^'^S) 




SMEC('^'^S) 


■^'^Sp'^S + ii] 


'^^Sp'T + ]-)] 


3.2909 


2.641 


3.690 


3.024 


3.822 



corrections for 0^*" and 2f states in "^^Ca and ^^S arc shown for selected daughter states. For 

01 state, the dominant energy contributions from couplings to the g.s. 1/2^ and the second 
excited state 5/2^ of A — 1 nucleus are almost equal ^. The dominant component of SMEC 
g.s. wave function is: 

[(4/2)'(«l/2)'] + [[(4/2)'(Sl/2)^]l/2+ {S)U] + (3) 
[[(4/2)'(M']5/2+]®(rf);c)]+--- 

, where index c denotes continuum state. Important contributions to E^^°'^'^^ are spread over 
six states in daughter nuclei. Contributions from couplings to the two lowest states {l/2f 
and 3/25*") favor ^^Ca state, whereas all contributions from couplings to higher excited states 
(7/2^", 3/2^, 5/2^ and 9/2^) favor ^^S state. As a result, the 2f state in ^^Ca is shifted 

2 The balance of E^T''^ in ^^Ca and ^^S is a direct consequence of the separation energies in these nuclei: 
5p(36Ca) = 2.56 MeV, Sn{^^S) = 9.89 MeV. 
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FIG. 4: The states Jf' in channel wave-functions [J^^' 



iJ"'^ = ^tX) which 



■"1 ' ^1 
36no /'SSr 



contribute most to the continuum-couphng correction from neutron continuum in ■^°Ca ( Ca states 
- red narrow boxes) and proton continuum in ^^S (^^P states - green wide boxes). The upper (lower) 
panel is for O^j^ {2'^) state. 
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down with respect to ^^S. The dominant component of the SMEC 2^ wave function is: 

m„2)\d,„f]y,. ® + (4) 

W„2)\s,/,)\d,/,)%,,^ ® {d)l^] + [[{d,/,f{sy,)\d,/,Y]y,. ® {d)l^] + 
W„2)\sy,)\„. ® (S);,)] + [[{d„2)\sy,)\d,/,)%,, ® {S)l^] +■■■ 

If one looks to the 2'^ — Q'^ energy difference in ^^Ca and '^^S, one may notice that a 
contribution to the energy of 0]*" from a couphng to l/2j'" in A — 1 nuclei is almost exactly 
compensated by contributions to the energy of 21 from couplings to 1/21/' and 3/2]*' states. 
The asymmetry in the position of 2^ states in '^^Ca and ^^S comes from the balance between 
the contribution of the 5/2^*' in the g.s. and contributions of higher lying states 7/2^, 
3/2^, 5/2^ and 9/2^^ in the 2^ state. These two sets of couplings act differently in the 
g.s. and in the first excited state. It is interesting to notice that the effect is produced by 
continuum couplings which practically do not involve s-wave neutron/proton components. 
These components, which are present both in the Oj'' state and in the 2^ state (cf Eq. (4)), 
mutually cancel out and do not contribute to the Thomas- Ehrmann shift in ^^Ca-'^^S mirror 
systems. 

The proximity of fp shell leads to important excitations from sd to fp across the = 20, 
Z = 20 shell closure. As a consequence, many new channels [J^'^~^;{l,i)Y^'''- become 
involved in the continuum coupling, leading to an enhanced spreading of contributions to 
the continuum-coupling energy correction. In general, an increased spreading of couplings 
leads to an increased fractionation of continuum-coupling contributions and the reduction 
of E'^'^°"^ due to an enhanced interference of large number of channels. 

Full SMEC calculations in sdfp shells for ^^Ca and ^^S are beyond actual possibilities so we 
)roceed by introducing the quenching factor Qf in the continuum-coupling energy correction 



5|]. SM calculations using WBT jl2|] and lOKIN [131] effective Hamiltonians, which include 
fp shell and allow for 2hjjj excitations, revealed that an appropriate quenching factor should 
be used to account for admixture of intruder configurations (cf Fig. [11]) • Almost identical 
admixtures of intruder configurations have been found both in parent nuclei and in all 
important states oi A — 1 nuclei. 

The value of the quenching factor depends on different extensions of S(i-shell effective 
interaction into a larger model space (cf Table 2). It turns out that the excitation energies of 
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TABLE II: Fraction of Ofko configurations in 2hL0 SM calculations for ^^Ca/^^S {A = 36) and 
corresponding daughter nuclei {A = 35). Results correspond to Of, 2^ states of the parent nucleus 
{A = 36) and to averages over 7 (resp. 8) most important states in = 5/2 (resp. = 3/2) 
daughter nuclei, respectively. 



Hamilt. 


space 


of {A = 36) 


2f{A = 36) 


{T. = 5/2),, 


{T. = 3/2),, 


lOKIN 


sdfp 


0.769 


0.768 


0.772 


0.782 


WBT 


sdfp 


0.787 


0.789 


0.795 


0.796 


WET 


psdfp 


0.687 


0.684 


0.683 


0.685 




0.5 0.6 0.7 0.8 0.9 1 

quenching factor 



FIG. 5: SMEC excitation energies of 2f states in ^"Ca and S are plotted as a function of the 
quenching factor. Experimental excitation energies of 2f state in ^^S and ^^Ca are depicted as 
diamond and circle, respectively. 

2f states in ^^Ca and ^^S are sensitive to the value of the quenching factor (see Fig. |5]). Both 
absolute excitation energies and their difference diminish as the quenching factor decreases. 
For Qf = 0.651, the SMEC excitation energy of 2f state reproduces exactly the experimental 
value in well-bound ^^S. This value for the quenching factor is close to the value found for 
WBT Hamiltonian (cf Table [Tll) in psdfp shells. 
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For lOKIN and WBT Hamiltonians in sdfp shells, one obtains higher values for Qf. For 
those values, SMEC reproduces well an experimental difference of 2f energies in ^^Ca and 
■^^S nuclei, but the absolute values of 2f energies are too high. The latter problem could be 
easily resolved by mirror-symmetry conserving correction of the Vq^^^-i^^ ^ monopole. 



V. CONCLUSIONS 



In this paper, we attempted to address two problems pertaining to the physics with 
exotic, weakly-bound nuclei: (i) What are the generic features of the continuum coupling in 
binding systematics of neutron-rich nuclei? (ii) Can one explain mirror-symmetry breaking 
effects in spectra by the asymmetry of threshold energies in mirror nuclei? 

The answers to these questions are not simple and require further investigations. Stud- 
ies in long chains of oxygen and fluorine isotopes revealed that the np residual coupling to 
the scattering continuum becomes strongly reduced with respect to the nn coupling in the 
vicinity of the neutron drip line. The np coupling is essential for an understanding of the 
anti-OES effect which is seen in odd-Z fluorine chain and leads to an apparent reduction of 
the gap parameter for neutrons. Hence, in weakly-bound nuclei close to the drip lines, the 
OES has three components: the first one originates from nucleonic pairing, the second one 
is the deformed mean-field effect jl^, and the third one originates from the np coupling via 
scattering continuum. One should stress that the latter component is of a completely differ- 
ent nature than the singular behavior of binding energies originating from the np correlations 



near N = Z line 



15|- 



Asymmetry in the spectra of mirror systems is another playground for the continuum shell 
model. In the mirror couple ^^Ca-'^^S, a relative shift of 2f excitation energies can be to a 
large extent directly related to the effect of the continuum coupling. However, in contrast to 
neutron-rich nuclei, the low one-proton separation energy in ^^Ca has little influence on the 
difference of excitation energies of 2^ states in ^^Ca and ^^S. This is due to the Coulomb 
barrier which suppresses continuum-coupling effects in weakly-bound systems close to the 
proton drip line. In that respect, weakly-bound/unbound systems at the proton drip line 
are radically different from those at the neutron drip line. To understand dissimilarity of 
nuclear systems at the proton and neutron drip lines is a challenge for the nuclear structure 
theory. 
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